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ABSTRACT 

We examine the fate of a dead radio source in which jet injection from the central engine has stopped 
at an early stage of its evolution {t = tj < 10® yr). To this aim, we theoretically evaluate the evolution 
of the emission from both the lobe and the shell, which are composed of shocked jet matter and a 
shocked ambient medium (i.e., shell), respectively. Based on a simple dynamical model of expanding 
lobe and shell, we clarify how the broadband spectrum of each component evolves before and after the 
cessation of the jet activity. It is shown that the spectrum is strongly dominated by the lobe emission 
while the jet is active (t < fj). On the other hand, once the jet activity has ceased {t > tj), the 
lobe emission fades out rapidly, since fresh electrons are no longer supplied from the jet. Meanwhile, 
shell emission only shows a gradual decrease, since fresh electrons are continuously supplied from the 
bow shock that is propagating into the ambient medium. As a result, overall emission from the shell 
overwhelms that from the lobe at wide range of frequencies from radio up to gamma-ray soon after 
the jet activity has ceased. Our result predicts a new class of dead radio sources that are dominated 
by shell emission. We suggest that the emission from the shell can be probed in particular at a radio 
wavelengths with the Square Kilometer Array (SKA) phase 1. 

Subject headings: particle acceleration — radiation mechanisms: non-thermal — galaxies: active — 
galaxies: jets — 


1. INTRODUCTION 

Relativistic jets in radio-loud active galactic nuclei 
(AGNs) dissipate their kinetic energy via interactions 
with surrounding interstellar medium (ISM) or intraclus¬ 
ter medium (ICM), and inflate a bubble composed of de- 
celer ated jet matter, which is often referred to as cocoon 
(e.g., Begelman et al.|I984 for review). Since the cocoon 
is highly overpressured a gainst the ambient ISM/ICM 
(Begelman & Ciofh 1989), bow shock is driven into the 
ambient matter. As a result, the shocked ambient gas 
forms a thin shell around the cocoon and envelopes the 
whole system. The thin shell structure persists until the 
cocoon pressure decre ases and the pressure equilibrium 
is eventually achieved ( Reynolds et al.pOOl l. 

The dissipation of jet energy due to the interaction 
with the ambient matter is accompanied by particle ac¬ 
celeration and non-thermal electrons are supplied to the 
cocoon from the jet. It is well established that prominent 
radio emission is produced by the accelerated electrons 
within the cocoon via synchrotron radiation. Hence, the 
cocoon is often termed as radio lobe. Inverse Comp¬ 
ton (IC) emission by the same population of electrons is 
sometimes observed at higher frequencies such as X-ray 


(e.g., Hardcastle et al.|2002| Isobe et al. 12002 

Croston et 

al. 2005 Kataoka &: Stawarz 120051) and, in some cases. 

even at gamma-rays (Abdo et al. 2010 

). 
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site of particle acceleration (iFuiita et al. 

2007 Berezhko 

2008 

Bordas et al. 2011 

Itoetal.|2011'K 

ino et al.|2013 ). 


rare and reported only at X-ray energies 

Kraft et al. 

2003 

Croston et al. 

|2007 

Jetha et al. 

2008 

. iVIoreoyer, 

most of them are tl 

rermal origin, anc 

only 

one source, 


the most nearby radio galaxy Centaurus A, shows evi¬ 


de nce for the non-thermal emission at X-ray (Croston et 


al. 2009). Due to the lack of observations, the number 


of studies that focus on the shell emission is small and 
therefore, its nature is poorly known. 

To tackle this issue, we have explored the evolution of 
the emission from both the lobe and shell simultaneously 
in the previous studies ( |Ito et aH201I Kino et al.|2013 ). 


Under the assumption that the jet is active throughout 
the evolution, it is shown that the overall spectrum is 
strongly overwhelmed by the lobe emission at most of 
the frequencies. Hence, detection of the shell emission is 
likely to be hampered by the lobe emission and limited 
in a narrow range of frequencies. This provides a natural 
explanation on the reason why the detection of shell is 
rare. 

While the large contrast between the lobe and shell 
emission is expected to be maintained as long as the 
jet continuously injects energy into the lobe, the situ¬ 
ation may change drastically after the jet activity has 


ceased. It is known from the literature (e.g.. 

Komissaroy 

& Gubanoy 1994 

Nath 2010 

Mocz et al. 2011 

[]) that lumi- 

nosity of the lobe 

emission s 

tiows a rapid decrease when 


the jet activity is stopped. This is because the radiative 
cooling leads to the depletion of non-thermal electrons, 
since fresh electrons are no longer supplied from the jet. 
Therefore, it is often claimed that these dead sources 
simply fade out soon after the jet activity ceases. How- 
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Ito et al. 


ever, the above statement may not be true when contri¬ 
bution from the shell emission is considered. Since the 
bow shock propagating into the ambient medium is ac- 


supplv fresh electrons 

in the shell (Reynolds & Begelman 

1997 

Reynolds et al. 

2002 

Ferucho et al. 2011). As a 


from the previous studies. 

First, let us summarize the basic assumptions in our 
model. (1) Regarding the ambient mass density, pa, a 
power-law dependence on radius, r, with an index, a, is 
assumed, Pa{?') = Here is the reference 

radius and po is the mass density at the radius. In the 


result, emission from the shell can persist its luminosity 
and eventually dominates over that from the lobe. How- 

present study, we adopt po = O.lrup cm ®,ro = 1 kpc 
and a = 1.5 based on the values inferred from the obser- 

ever, up to present, there are no quantitative studies that 

vation of elliptical galaxies (e.g.. 

Mulchaey & Zabludoff 

focus on this issue. 

19981 IMathews & Brighenti|2003 

h'ukazawa, iVlakishima, 

Regarding the evolutionarily track of radio sources, it 

& Ohashi 

2004 [lukazawa et al. 

2006|). Here rUp is the 


i.e.. Gigahertz Peaked Spectrum (GPS) and Compact 
Steep Spectrum (CSS) so urces, represen t the young stage 
(- 10^ - 10® yr; e.g.. 


|1998) of typical well- 
lat have estimated ages 


O’Dea 

extended classical radio galaxies t 

of - 10^ - 10® yr. On the other hand, it is not clear 
whether “all” compact radio sources can survive to the 
ages of typical radio galaxies. This comes from observa¬ 
tional fact that these compact sources represent a large 
fraction (~ 15 — 30%) in the flux-l imited catalogue of the 


radios sources (Fanti et al. 1995), which is much larger 
than the expected value (^ 0.1%) from their youthness. 
A simple and plausible explanation for this discrepancy is 
that significant fraction of young compact radio sources 


(5, 10^ yr) ( Alexander||2000 Marecki et al.||2003 Gugli- 

ucci et al. 


Urienti & Dallacasa||2U10| Urienti et 

al. 

2010 Kunert- 

lajraszewska et al. 2010). Then, large 


population of compact dead radio sources in which jet 
injection from central engine has ceased is expected to 
be hidden in the universe. Therefore, clarifying the evo¬ 
lution of dead radio sources including the contribution 
of shell is essential for revealing the evolution of radio 
sources. 

Motivated by these backgrounds, we explore the emis¬ 
sion from radio sources in which jet activity has ceased 
at early stage of their evolution in the present study. In 
particular, we focus on the evolution of the relative con¬ 
tribution of the lobe and shell emission. We show that 
the shell will be dominant at most of the frequencies from 
radio up to gamma-ray soon after the jet is switched off 
and discuss the possibility for detecting the emission. 

This paper is organized as follows. In Section!^ we in¬ 
troduce our dynamical model, which describes the evolu¬ 
tion of the shell and cocoon, and explain how the energy 
distribution of the electrons residing in these regions and 
the spectra of the radiation they produce are evaluated. 
The obtained results are presented in Section The 
summary and discussion of our results is given in Sec¬ 
tion 1^ Throughout the paper, a AGDM cosmology with 
Hq = 71 km s~^ Mpc~^, Hm = 0.27 and Ha = 0.73 is 


adopted (Komatsu et al. 2009). 


2. MODEL 


Following the previous studies (Ito et al. 2011 Kino 


eTaTjlM^ , we evaluate non-thernial emissions from co¬ 


coon and shell based on a simple analytical model that 
describes their dynamical expansion. The model is iden¬ 
tical to that employed in the previous studies during the 
early phase when the jet is active. The only difference 
is that we properly take into account the evolution after 
the jet injection has ceased. In this section, we briefly 
review the employed model and explain the differences 


p- 


proton mass and po.i = Po/O.im' 
coon and shell are approximated as a sphere and neglect 
its elongation along the jet direction merely for simplic- 
ity. (2) We use the stan dard thin shell approximation 


(Ostriker & McKee 1988) and assume that most of the 
matter swept up by the bow shock is concentrated in 
a region of width , 5R, behind the shock which is thin 
compared with the radius of the shock, R. 

Regarding the dynamics of expansion, we consider two 
phases depending on the age, t: (i) the early phase in 
which the jet injection is present (t < tj) and (ii) the late 
phase in which the jet injection has ceased {t > tj). Here 
tj denotes the duration of the jet injection. As mentioned 
abo ve, we use the same model used in the previous stud¬ 
ies (Ito et al. 2011 Kino et ^|2013) in the early phase 


(i? = dR/dt oc After the energy injection 

from the jet ceases, the cocoon will rapidly lose its en¬ 
ergy due to adiabatic expansion and give away most of 
its energy into the shell within in a dynamical timescale 
^ R/R. Hence, after the transition time ~ R/R ^ tj, 
cocoon pressure becomes dynamically unimportant, and 
the energy of the shell becomes dominant. Therefore, in 
the late phase, expansion of the bow shock is expected 
to akin to Sedov-Taylor expansion {R oc t“(3-Q;)/(5-a)^_ 
Based on these considerations, we can approximately de¬ 
scribe the two phases continuously by connecting the the 
expansion velocity of the bow shock at t = tj 


as 


. -(2-a)/(5-a) 


Ro (t') for 0 < t < tj, 

= <!. ; -( 3 -«)/( 5 -«) (1) 

Ro (^jrj for tj < t. 


where 


Ro = C^ 


\PorS 


l/(5-c 


^(2-a)/(5-Q) 


( 2 ) 


Here Lj is the power of the jet and C = 3/(5 — 
a) [(3 — a) (5 — a)®(// — l)/{47r(2a^ -I- a — 18//a -I- 63// — 
28)}]^/^®““\ where 7 c is the specific heat ratio of the 
plasma inside the cocoon. From the above equation, 
the radius of the bow shock is determined as R{t) = 
J* R{t')dt'. In the present study, we assume % = 4/3, 
since plasma within the cocoon is expected to be rela¬ 
tivistic. 

Once the expansion velocity R is obtained, the prop¬ 
erties of the shell are determine d in the same m anner 
described in the previous study (Ito et al. 2011). We 


employ the non-relativisti c Rankihe-Hugoniot cond ition 
in the strong shock limit (Landau & Lifshitz 1959) and 
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evaluate the density and pressure in the shell as 

Ps{t) = ^ ] ps.{R), 

7a- 1 


and 


Psit) = 

7a + 1 


(3) 


(4) 


respectively, where % is the specific heat ratio of the 
ambient medium. In the present study, we adopt % = 
5/3 since the temperature of the ambient medium is non- 
relativistic. The shell width, SR, in which most of the 
mass swept up by the bow shock is contained is evaluated 
from the conservation of mass psVs = /q pa,{r)dr, 
where 14 = ^ttR^SR is the volume of the shell. Hence, 
the shell width is given by SR = (ya — l)[(7a + 1)(3 — 
As a result, the total internal energy of the 
shell, Es = Psl4/(7a — l)j scales linearly with time as 
Es = fiLjt oc t for t < tj and gradually asymptotes 
to a constant value of Eg = f 2 Ljtj oc for t tj. 
Such a behaviour is naturally expected, because energy 
is continuously injected into the system from the jet with 
a constant rate when the jet is active, while there is no 
energy source after the jet activity ceases. Here /i and 
/2 are the fractions of energy injected by the jet that are 
converted into the internal energy of the shell which are 
given by /i = 18(7c - 1)(5 - a){% + l)~‘^[2a^ + (1 - 
187c)q; -I- 637c — 28]“^ and /2 = (3/2)^““/i. For the 
values employed in the present study, ( 7 c = 4/3, 7 a = 
5/3 and a = 1.5), /i ~ 0.11 and /2 ~ 0.21 are obtained. 
Hence, the ratio of the internal energy of the shell to the 
total energy deposited in the system is roughly constant 
10 — 20%) throughout the evolution. 


As described in the previous study (see § 2.1 of Ito et al. 


2011 


for detail), for t <tj, the pressure of the cocoon is 
directly obtained from the dynamical model and is given 

by 


Pc(t) = fpP,{t), 


(5) 


where fp is the ratio of the pressure of the cocoon to that 
of shell and is given as fp = ( 7 a -f 1)(7 — 2a)/[6(3 — a)]. 
For typical numbers, 7 a = 5/3, 0 < a < 2, fp depends 
on a only weakly and /p ~ 1. This is expected, since 
the pressure o f the shell and co coon should be roughly 
equal (see e.g., Heinz et al. ]I9M1 ). Since pressure balance 


is also expected for t > tj, we assume that the ratio of 
the pressure is maintained and determine Pc from the 
same equation also in the later phase. Regarding the 
evolution of the radius of the cocoon, Rc, we assume 
Rc = R ior t < tj as in the previous studies. On the 
other hand, since the cocoon expands adiabatically, the 
cocoon radius satisfies 744/7“ = const, for t > tj, where 
Vc = 47 ri ?^/3 is the volume of the cocoon. Hence, the 
evolution of the cocoon radius can be summarized as 

' R{t) for 0 <t <tj, 

, 1/37c . (6) 


Rcit) = 




As a result, cocoon expands slower than the shell in 
this phase and makes the region of the shocked shell 
wi der. It is worth noting that sim ilar behaviour is found 
in Reynolds & Begelman (1997). The total internal 


energy deposited in the cocoon is evaluated from the 
obtained pressure and radius as Ec = Pc 44 /(7c — !)• 
For t < tj, Ec scales linearly with time as Ec = 
f\ohe,iLjt as in the case of the shell, where /iobe,i = 
(5 — a)(7 — 2a)[2Qf^ -|- (1 — 187c)q! + 637c — 28]“^. At 
later time {t > tj), the evolution of Ec asymptotes to 
Ec = fiohe, 2 Ljtj{t/tj)~^^'^-~'^P^'^^^^~°‘'>\ where /iobe .2 = 
(2/3)“(7c-i)/7cP or the values employed in the 
present study, (% = 4/3, 7 a = 5/3 and a = 1.5), 
fiohei 0.54 and /iobe ,2 ~ 0.46 are obtained. 

In determining energy distribution of non-thermal elec¬ 
trons, N{'^c,t), within the shell and lobe, we approxi¬ 
mate each region as one-zone and solve the kinetic equa¬ 
tion which takes into account the injection and cooling 
of particles given by 

= ^[7cool(7e,i)A^(7e,0] +Q(7e,i). (7) 

Here 7 e, 7cooi = —d'je/dt, and Q are the Lorentz factor, 
the cooling rate via adiabatic and radiative losses, and 
the injection rate of accelerated electrons, respectively. 
The injection rate and the cooling rate are determined 
from the dynamical model in the same manner as de¬ 


scr ibed i n the previous studies (|lto et al.||2011 
Regardms 


Kino et 


al 


.egarding the shell, we assume a continuous injection 
throughout the evolutions given by 

Qile,t) = K{t)-ffP for O^min Te ^ Tmax it). ( 8 ) 

Here 7 min and 7 max are the minimum and maximum 
Lorentz factors of the accelerated electrons, respectively. 
In the present study, we employ ymin = 1 and p = 2. 
In determining the maximum Lorentz factor, we assume 
that the acceleration is limited by cooling. As for the 
electron acceleration mechanism, we assume the well- 
known diffusive shock acceleration in which the acceler¬ 
ation rat£_ca£_^ written as 7accei = ieBJP/{20fnieC^) 
(e.g.,|Drury|1983|). Here B is the magnetic field strength 
in the shell and ^ is the so-called “gyro-factor” which can 
be expressed as the ratio of the energy in ordered mag¬ 
netic fields to that in turbulent ones (^ = 1 corresponds 
the Bohm limit). Hence, the maximum Lorentz factor 
is given by the value where cooling rate, 7 cooi( 7 e), and 
acceleration rate 7 accei become equal. The normalisation 
factor, K, is computed by assuming that a fraction, Cg, of 
the shock-dissipated energy is carried by the non-thermal 
electrons. In our calculation, we assume constant energy 
injection rate throughout the evolution, and, for a given 
source age oit = Ggej the normalization factor at an arbi¬ 
trary time {t (< tage))) is determined from the equation: 

= eePs(tage)/tage- We 

choose C = 1 and and = 0.01 as fiducial values for the 
parameters that characterize the acceleration efficiencies 
based on the ob s ervations of supernova remnants (e.g., 
Dyer et allpOOl |Ellison et al.||2001 |Bamba et al.||2003[ 


Yamazaki et al. |2004| |Stage et al.| 2006| |Tanaka et al.| 

^008| since the properties of the shock are similar to 
that of the bow shock considered here. 

^ We note that the tempor al evolutions of and Q are taken 
into account in this paper an diKino et al.|||2013ll. whereas they were 
neglected in |Ito et al. ]poII|i. -- 
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Considering the electron injection into the lobe, we as¬ 
sume continuous injection as in the case of the shell for 
t < ty In contrast, it is assumed that injection of the 
non-thermal electrons ceases for t > ty since the jet can 
no longer supply matter in this phase. Hence, the en¬ 
ergy distribution of the electrons injected into the lobe 
is given as 

r (0 < t < tj) 

Qlobe('Te;^) — N Tmin,lobe <7e< O^maXjlobe ? (9) 

lo 

where 7min,iobe and 7max,iobe are the minimum and max¬ 
imum Lorentz factors, respectively. The imposed values 
for the power-law index and minimum Lorentz factor are 
the same as those adopted for the shell (piobe = 2 and 
7min,iobe = !)■ On the other hand, we use as fixed value 
7 max,iobe = 10^ for the maximum Lorentz factor. The 
employed value of 7 max,iobe is based on multi-wavelength 
observations of FRII radio galaxies which suggest elec¬ 
trons within the lobe do not have Lorentz factor well be- 


yond ~ 10 

(IStawarz et al. 

2007 

Yaji et al. 

2010p. As shown in 


Godfrey et al. 
noted, 


it IS 

ever, tJrat tne conclusion ot tJre present study is in¬ 
sensitive to the assumed values of piobe, 7min,iobe and 
7max,iobe- In determining itiiobe, we assume that a frac¬ 
tion, Cejobe, of the energy deposited in the lobe is carried 
by non-thermal electrons and the energy injection rate 
of the non-thermal electrons are constant up to t = tj: 
/7""(7e - l)™ec2Q(7e,brf7e = ee.lobe^^c (tj) Aj = 
/lobe.iTj. Although the value of Cejobe is not constrained 
very well, we adopt Ce.iobe = 1 in the prese nt study which 
is often assumed in previous studies (e.g. , iStawarz et al. 


2009^ 

howl 


2008 Ostorero et al. 


1[m0l[!toiraI1[MT|. 

loss rate, 7 rnoi, we take 


Regarding energy loss rate, 7coob we take into account 
radiative and adiabatic (expansion) losses. The adiabatic 
cooling is evaluated from the expansion of the system as 
7 ad = {R/R)le for the shell and 7ad = {Rc/Rc)le for the 
lobe, where i?c = dRc/dt. As for the radiative cooling, 
we consider the synchrotron radiation, 7 syn, and the IC 
scattering 71 c. 

The typical magnetic field strengths in elliptical galax¬ 
ies and c lusters of galaxies are inferred to be around few 
liG (e.g.. Moss fc Shukurov|1996 [ Vikhlinin, Markevitch, 


fc MurravpOTT 

Schekochinin et 


Clarke et ai. 


al J2005p . 'T 


S -GarlllifcTM 

ire, we adopt magnetic 


TIeId~sfrengnr'br"B''^^T 0 /rG in the shell for the calcu¬ 
lation of synchrotron emission, since shock compression 
leads to amplification by a factor of ^ 1 — 4. On the 
other hand, magnetic field strength within the lobe. Be, 
is determined under the assumption that a fraction eb, 
of the energy is carried by the magnetic fields. In 
the present study we assume eb = 0 . 1 , a magnetic field 
strength factor of a few below the equipartition value 


(e.g., Isobe et al. 2002 Kataoka et al. 2003 

Croston et 

al. 20 (W 

2005 Kataoka & Stawarz||2005). 


Regarc 

ling the 1C scattering, lull Klein-IN ishina (KJN) 


cross section is taken into account (Blumenthal & Gould 
1970P and, as a source of seed photon helds, we con- 


and cosmic microwave background (CMB). We modelled 
the spectra of the photons from the disc, torus and host 
galaxy with a black-body spectra peaking at frequencies 
given by v\jy = 2.4 x 10^^ Hz, = 1-0 x 10^^ Hz and 
t'NiR = 1.0 X 10^"^ Hz, respectively. As for the luminosities 
of the emissions, we adopt Lyv = L tr = 10 "^^ erg s~ ^ 


for the disc and torus emission (e.g., Elvis et al. 1994 
Jiang et al.||2006|, and Lntr = 10'^^ erg s~ f or tfe host 
galaxy emission (e.g., de Ruiter et al. 20051. GMB is 
given by a black-body that has peak frequency and en¬ 
ergy density given by j^cmb = 1-6 x 10^^(1 -|- z) Hz and 
Gcmb ~ 4.2x 10 “^^(l-|- 2 :)^ erg cm“^, respectively, where 
z is the cosmological redshift of the source. The spectrum 
of the seed photon field of the synchrotron emission from 
the lobe is obtained in our calculation self-consistently. 
The IC cooling is determined by the given spectra and 
source size as in the previ ous study (for detail, see § 2.2 


and §4.1 of Ito et al. 2011 1 . 

The energy distribution of non-thermal electrons, 
A^( 7 e,t), is obtained by putting the evaluated injection 
rate, Q( 7 e, t), and the cooling rate, 7 cooi( 7 eji), in Equa¬ 
tion Q . The spectra of the synchrotron and IC emission 
are s^-consistently calculated from the obtai ned energy 


distri bution (for detail, see §2.3 and §4.1 of Ito et all] 


20111 . 

As shown in the previous studies, thermal emission is 


emission (e.g., Heinz et al. 

1998 Kaiser & Alexander! 

1999| 

Reynolds et al.||2001| |Z,anni et al.||2003t |Bordas| 

et al. 

201ip. Therefore, we 

also calculate the thermal 


study. Under the assumption that most of the shock 
dissipated energy is converted into that of the thermal 
electrons, the luminosity is estimated as L,/,brem = 
2^TTe^ 

where is the Boltzmann constant, and rig = Psjmp, 
and T = Ps/2nefcB are the number density and tem¬ 
perature of the thermal electrons within t he shell, 
respectively. Here gff is the Gaunt factor (jRybidd 
& Lightman 1979) and is set to be unity merely tor 
simplicity. 

3. RESULTS 

In this section, we show the time evolutions of the en¬ 
ergy distribution of non-thermal electrons and the re¬ 
sulting emission. As a fiducial case we focus on sources 
with jet power of Lj = 10^^ erg s“^. Since our aim is to 
explore the evolution of short-lived sources in which the 
energy injection into the lobe have stopped in the early 
stage of their evolution, we set the duration of energy 
injection as tj = 10® yr. The chosen value of the age is 
around the upper end o f the estimated age of the com¬ 
pact radio source (e.g., Murgia 2003[ |Polatidis fc Con- 


way 


20031). The overall linear size ai t = t 
0.86 / 0 (C A' where L 45 = Lj/lO'^® erg s 

tj. 5 =tj/ 10 ® 


and 


yr- 


It is noted that, since spherical symmetry is assumed, 
the source size and age correspondence predicted in our 
model show relatively large discr epancy from those in- 


sider U V emission from the accretion disc, IR emission 

Taylor et al.||2000 

Conway 

2002 ). 

from the dusty torus, stellar emission in NIR from the 

ometry is elongated along 

the jet 


host galaxy, synchrotron emission from the radio lobe 


1999 

I'he actual source ge- 
he jet axis, and, for a given 
source age, their linear sizes along the jet direction and 
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transverse direction are larger and smaller by a factor 
^ a few, respectively. It is emphasized, however, that the 
discrepancy does not have signihcant effect on the essen¬ 
tial features of our results, since large axial ratio is rarely 
observed. Moreover, the axial ratio is reduced after the 
jet activity ceases {t > tj), i.e. the phase that we mainly 
focus on this study, since the thrust from the jet which 
causes the elongation is no longer present. Therefore, 
we can consider that the results obtained in the present 
study at a given age t apply to the sources whose linear 
extension is larger than R{t) within a factor of few. 

3.1. Evolution of Non-thermal Electrons 

In Figure we display the energy distribution of non- 
thermal electrons within the lobe {left panel) and shell 
{right panel) for different source ages {t = 10^, 1.1 x 10®, 
2 X 10®, 5 X 10® and 10® yr). When the jet is active {t < 
tj), the non-thermal electrons are continuously supplied 
in the shell and lobe. In this initial phase, the energy 
distributions of electrons within the lobe and shell show 
a broken power-law shape. The break is located at the 
energy where the cooling time of the electrons becomes 
roughly equal to the dynamical time {t ~ GooKbe))- 
Below the break energy, the electron energy distribu¬ 
tion maintains the form of the injected energy spectrum 
which is roughly given by N{-je) ~ Q{le)t oc since 
the cooling effect is negligible at these energies. Above 
the break energy, the energy distribution is modihed by 
the cooling loss and can be roughly approximated as 
N{'ye) Q(7e)Gooi(7e) with a sharp cut-off at the in¬ 
jected maximum energy (ymax and 7max,iobe, for the shell 
and lobe, respectively). While for the shell the dominant 
radiative cooling is provided by the IC scattering off the 
disc and torus photons, synchrotron cooling is slightly 
larger than the IC cooling for the lobe. Since the cooling 
is more rapid in the l obe, the break energy is lower than 
that of the shell (see Ito et al.||2011 for more detail). 

After the jet activity has ceased {t > tj), the evolu¬ 
tion of electron energy distribution differs largely be¬ 
tween the shell and lobe. Regarding the shell, the evolu¬ 
tion does not show large change from that at the early 
phase {t < tj). The electron energy distribution can be 
roughly described by a broken power-law shape. The ra¬ 
diative loss is mainly dominated by the IC scattering of 
the torus photons. Since the energy density of the pho¬ 
ton decreases with the source size Uph oc R~^, the break 
energy gradually increases with the age. 

On the other hand, due to the cessation of the elec¬ 
tron injection {Q{'je) = 0), the evolution of the energy 
distribution of the electrons in the lobe changes drasti¬ 
cally from that at the early phase {t < tj). Since fresh 
electrons are no longer supplied, population of electrons 
with energy above the break energy determined at < = tj 
rapidly decreases because of the cooling loss. The cool¬ 
ing proceeds predominantly through synchrotron emis¬ 
sion and the corresponding break energy is found at 

7br,lobe ~ ll^eC/(4(T7’17B,lobeO 

^4.8 X (10) 

where t/B.iobe = -i = cb/O.I and = 

tj/10®yr. The rapid depletion continues until t ~ 2tj at 
which all the high energy electrons (7e > 7br,iobe) are 


cooled. Thereafter, electrons gradually cool mainly due 
to adiabatic loss ( 7 e oc oc t~^By These features are 
clearly seen in Figure 

3.2. Evolution of Broadband Emission 

In Figure we show the corresponding evolution of 
the total photon fluxes, rf^, from the shell (solid lines) 
and lobe (dashed lines) for a source located at the red- 
shift of z = 1. Together with the total photon fluxes, 
the contributions from the synchrotron emission, and 
IC scattering of various seed photons are displayed in 
Figure for a source with t = 10® yr and 10® yr. As 
mentioned in the previous section, contribution from the 
thermal bremsstrahlung emission is also considered for 
the shell emission. Attenuations due to synchrotron self¬ 
absorption (SSA) and pair production due to photon- 
photon interaction are also taken into account in the 
spectra. The absorption coefficient for SSA within the 
lobe, ckiobe,!/, and shell tts.y are determined independently 
from the magnetic held strength and el ectron energy dis¬ 
tribution w ithin these regions following [Rybicki fc Light^ 


man 


(1979). The spectra of the shell and lobe below the 


frequency which satisfy a^^^R = 1 and aiobe,;vRc = 1 (op¬ 
tically thick), respectively, are described as f^ cx 
Note that the absorption in the lobe is much stronger 
than that of the shell since the number of non-thermal 
electrons and magnetic held strength within the lobe is 
larger than those of the shell. Gamma-rays above energy 
> 100 GeV are subject to pair production due to photon- 
photon interaction with optical-infrared background ra- 
diation . Here we used the model of IFranceschini et al.l 
(2008) for the background radiation and ejraluat^ the 
attenuated photon huxes shown in Figures and We 
also plot the unattenuated gamma-ray spectra in Fig¬ 
ure [3 (dotted line). 

3.2.1. Lobe emission 

Regarding lobe emission, synchrotron and IC compo¬ 
nents produce prominent emission extending from radio 
up to hr ~ 10 GeV when the jet is active {t < tj). Their 
overall luminosities are roughly comparable and propor¬ 
tional to energy injection rate into the non-thermal elec¬ 
trons {uL^ oc Ce^ioheLj). This is simply because the en¬ 
ergy injected into the non-thermal electrons Ce^ioheLj is 
immediately converted into the radiation. 

After the jet activity ceases {t > tj), due to the absence 
of particle injection, radiative cooling leads to rapid de¬ 
pletion of high energy electrons 7 e > 7 br,iobe within the 
lobe. As mentioned earlier, this rapid decrease contin- 

radiate 
As a 


ues until all the electrons above the break energy 
away most of their energy t ~ 2tj (Figure, [l 
result, the high energy cut-off frequency rapidly shifts 
to the lower energies and the overall emission fades out 
quite rapidly, as seen in Figure It is worth noting 
that this passive feature of the l^e emissi on is consis- 


& Gubanov 

1994 

Nath 

o 

1—1 

o 

(N 

Mocz et al. 


2 Ull| and 


ing radio source s, so-called faders, are rarely found in 
the observations ( Kunert-Bajraszewska et al.|2005 2006 


Orienti et al.lpOlII ). After this passive phase {t N 2tj j , 
although relatively slow, the emission continues to fade 
predominantly through adiabatic cooling. It is stressed 
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Figure 1. Energy distributions of non-thermal electrons within the shell {right panel) and lobe {left panel) for source with jet power of 
Lj = 10'*®erg and injection duration of tj = 10®yr. The various lines display the cases for sources with ages of 10® (red line), 1.1 X 10® 
(blue line), 2 X 10® (purple line), 5 X 10® (light blue line) and 10® yr (black line). 


that the evolution at the fading phase (t > tj) is insensi¬ 
tive to the assumed value of 7max,iobe or the assumptions 
on the seed photon fields, since 7br,iobe does not depend 
on these assumptions. 

3.2.2. Shell emission 

Regarding emission from the shell, as in the case of the 
lobe, the overall luminosity of the non-thermal emission 
is proportional to the jet power and the energy fraction 
of the non-thermal electrons {vLi, oc eeLj) when the jet 
is active t < tj. However, much fainter emission is pro¬ 
duced in the shell, since smaller energy fraction is as¬ 
sumed (ce <SC Ce.iobe)- Unlike in the case of the lobe, IC 
emission is brighter than the synchrotron emission, since 
the energy density of magnetic field, Ub = is much 

smaller than that of the seed photons. (Ub ^ Uph). 
Since the maximum energy of the electron is larger than 
that of the lobe (7max ^ 7max,iobe), the overall spectrum 
extends up to higher frequencies. The high energy cut-off 
is observed at /i;/ ~ 1 — 10 TeV, and at this high energy 
range the spectrum is modified by the attenuation due 
to pair production ( Figur e^. 

As mentioned in §3.1[ tne evolution of non-thermal 
electrons at the later phase {t > tj) does not show large 
change from the initial phase {t < tj). The overall lumi¬ 
nosity continues to be dominated by the IC emission. As 
in the case of the lobe, the emission from the shell be¬ 
gins to fade in this phase. On the other hand, however, 
their cut-off frequency /iz/ ^ 1 — 10 TeV remains nearly 
unchanged and the decrease in the luminosity proceeds 
much more slowly (Figure]^. As in the initial phase, 
the overall luminosity is proportional to the energy injec¬ 
tion rate into the non-thermal electrons {i^L^ oc e^Ljtj/t). 
Hence, after the jet injection has ceased the luminosity 
gradually decreases with age as vLv oc if the energy 
fraction of electrons Ce doe not vary with time. This 
large difference in the behaviour of luminosity evolution 
between the lobe and shell is mainly due to the fact that 
the fresh electrons that are accelerated at the bow shock 
are continuously supplied into the shell even after the jet 
cessation. 

In contrast with the non-thermal emission, the ther¬ 
mal bremsstrahlung emission from the shell has a rel¬ 


atively weak dependence on the jet power and has a 
strong dependence on the ambient medium. The emis¬ 
sion produces a bump in the spectra at X-ray frequen¬ 
cies {hv ~ 1 — 100 keV) with peak frequency given as 
hv ^ UbT cx R(t)‘^ and dominates over the non-thermal 
emission as seen in the Figures. [^and[^ 


3.2.3. Comparison between lobe and shell 

Finally, let us compare the emission from the lobe and 
shell. In the initial phase when the jet is active (t < tj). 


the result does not var 

y from those obtained 

in the pre- 

vious studies ( 

Ito et al. 

2011 

Kino et al.|2013 

), since the 

same set up is used tor the calculation. 1 

he lobe emission 


dominates over the shell emission up to their cut-off fre¬ 
quency {hv ^ 10 GeV). At higher energies, however, the 
shell emission becomes dominant up to tiz/ ~ 1 — 10 TeV. 

On the other hand, soon after the cessation of the jet 
activity {t > tj), emission from the shell becomes dom¬ 
inant at most of the frequencies due to the rapid de¬ 
crease of the lobe luminosity. For example, the emission 
from the shell overwhelms that of the lobe at frequencies 
> 10^ Hz for sources with age of t > 10® yr (see Figure]^. 
As a result, dead radio sources eventually become domi¬ 
nated by the shell emission even at the radio frequencies. 
This shell-dominated phase is expected to continue up to 
the age where the expansion velocity becomes subsonic, 
since particle acceleration is unlikely to take place there¬ 
after. As discussed in next section, the corresponding 
age can be as large as > 10^ yr. Until then, luminosity 
of the shell gradually decreases as oc It is noted that 


vations (Kunert-Bajraszewska et al. 

200 ^ 

2006 

Orienti] 

et al. 

20TU 

) should be in passive p 

lase {i 

‘j < t 

<^tj) 


but still dominates the emission at radio wavelengths. 


4. SUMMARY AND DISCUSSIONS 

In the present study, we have explored the evolution 
of the emission from the lobe and shell of radio sources 
which have ceased their jet activity in early stage of their 
evolution {tj < 10® yr). It is shown that the lobe emis¬ 
sion rapidly fades after the jet injection has been stopped 
{t > tj) due to the absence of particle supply from the 
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t=10^yr 

t=l.lxl0^yr 

t=2xl0V 

t=5xl0^yr 

t=10®yr 


Figure 2. Broadband spectrum of source with a jet power of Lj = 10^®erg s“^ that located at the redshift of z = 1, and an injection 
time tj = 10® yr. The dashed lines and solid lines display the emission produced within the lobe and shell, respectively, for sources with 
age of 10® yr (red line), 1.1 X 10® yr (blue line), 2 X 10® yr (purple line), 5 X 10® yr (light blue line) and 10® yr (black). 


jet. On the other hand, shell emission only shows grad¬ 
ual decay even after the jet activity has stopped since 
supply of electrons from the bow shock is maintained. 
Therefore, although faint, emission from the dying ra¬ 
dio sources eventually becomes shell dominated at wide 
range of frequencies from radio up to gamma-ray. 

Our results suggest that shell emission is essential for 
studying the nature of dead radio sources. Shell emission 
may be detectable for exc eptionally high p ower sources 


(Lj > 10 ^® ergs e.g., Ito et al. 2008) that are lo¬ 


cated nearby {z < 0 . 2 ). On the other hand, for typical 
sources (Lj ~ 10^^ erg s“^) located at redshift of z ~ 1, 


a ty pical value for compact radio sources; e.g., |0’Dea| 


1998), the emission is below the thresholds of Rie cur¬ 
rent instruments at any frequencies (see Figure.]^. This 
provides a natural explanation why such a shell emission 
has not been reported so far. Regarding the future mis¬ 
sions, SKA telescope is capable of detecting the emission 
in the radio band. As shown in Figure.]^ the detection is 
marginal for SKA phase 1 (the emission slightly exceeds 
3cr detection limit for an integration time of 10 hours with 
conservative values of Ce = 0.01 and Lj = 10 '^® erg s“^). 
It is noted that the emission is roughly linearly propor¬ 
tional to the acceleration efficiency and jet power Lj. 
Since these values inevitably have dispersion from source 
to source, we can expect a certain fraction of sources to 
be well above the detection limit. Moreover, in the phase 
2 (SKA2), the sensitivity is expected improve by an or¬ 
der of magnitude. Hence, SKA will be a powerful tool 
to reveal the population of dead radio sources which are 
dominated by the shell emission. 

In order to confirm that the emission comes from the 
shell rather than lobe, spatial angular resolution of the 
observation is important. If the emission is indeed pro¬ 
duced in the shell, we expect to see limb brig htening fea¬ 
tures (|Heinz et al.||1998| Bordas et al.||2011|). The emis¬ 
sion will appear as an extended source with a total an¬ 
gular size given as ~ 2{2R/5 kpc)(I?A/I Gpc)“^ arcsec. 


where Da is the angular diameter {Da ~ 1-6 Gpc for 
z = I). Hence, high spacial resolution (~ 0.1 arcsec) 
provided by SKA can detect such a surface brightness 
distribution. 

The shell emission dominated phase is expected to 
persist until the expansion velocity R{t) becomes sub¬ 
sonic, since, thereafter, strong shock will no longer be 
present. The sound speed in ambient gas is given as 
Cs = [la.kBTa./{iJ,mp)Y^^, where Ta and /r = 0.6 are the 
temperature and mean molecular weight of the ambient 
gas, respectively. Then the age when the source becomes 
subsonic {R/cg = 1) can be estimated as tg « 8.4 x 

yr, "where '-L — feu 'I I, 1 ke\^. 
From the above equation, it is confirmed that the dura¬ 
tion of shell dominated phase is much longer than that 
of the lobe dominated phase {tg — Therefore, if 

large fraction of young radio sources die before becoming 
classical extended radio gala xies (t ^ 10 ^ — 10 ^ yr) as in- 
dicated in previous studies (|Alexander||2000l iMarecki et 

& Dali 


al.|2003[ Gugliucci et al.pOOSUGrienti & Dallacasal 


Orient! et al.|2010 Kunert-Bajraszewska et al.|2010p, our 

results predict the existence ot a large population of un¬ 

detected shell dominated sources residing in the universe. 
We expect that future survey by SKA will enable us to 
perform systematic study of these sources and quantify 
the fraction of dying sources through the detection of 
their shell emission. These observations will provide an 
important constraint on the evolution of radio sources. 

Lastly, let us comment on the radio sources which stop 
their j et activity after becomi ng classical extended radio 
lobes ( Fanaroff fc Riley||1974 ). Although we focused on 
the sources which die at and early stage of their evolution 
(< 10 ® yr), qualitatively similar results are expected for 
such sources. As shown in the present study, their lobe 
emissions fade rapidly and eventually become shell domi¬ 
nated. Actually, candidates of dead classical radio galax¬ 
ies that are in the lobe fadying phase have been found in 
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Figure 3. Broadband spectrum of source with a jet power of Lj = lO^^erg s ^ that is located at the redshift of 2 : = 1. The top and bottom 
panels display the spectra for sources with ages of 10^ and 10^ yr, respectively. The thick black and thick red lines show the total flux from 
the shell and lobe, respectively. The various thin dashed and thin solid lines show the contributions to the total emission produced in the lobe 
and shell, respectively, from the synchrotron emission (red line) and the IC scattering of UV disk photons (light blue line), IR torus photons 
(orange line), NIR host galaxy photons (blue line), CMB photons (green line) and lobe synchrotron photons (purple line). The contribution 
from the thermal bremsstrahlung emission is also shown for the shell (yellow line). The dotted lines show the case when the attenuation 
due to pair production is neglected. Also shown are the sensitivities of the XMM-Newton, ALMA, Jansky VLA and SKAl (SKA1_MID 
and SKAl_LOW). The assumed integration times for XMM-Newton, Jansky VLA and ALMA are 100 ks, 4 and 10 hours, respectively. The 
sensitivity of ALMA is calculated by using the ALMA Sensitivity Calculator (https://almascience.eso.org/proposing/sensitivit v-calculator). 
As for SKA, we display 3a detection limit for an integration time of 10 hours which is calculated based on the Appendix of|Prandoni &| 
|Seymou7| | |2014[ |. 


the observations (e.g., Murgia et al.|20lT Hurley-Walker 


et al.|[2M5|. Therefore, we also expect to hnd the scaled 


up version of the compact shell dominated sources con¬ 
sidered in the present study. Such sources are also very 
interesting and detectable by the SKA. Provided that 
the observed numbers of compact {t < 10 ^ yr) and well 
extended radio so urce {t ^ 10 ^ — 10 ® yr) are comparable 
(|Fanti et al.|1995|) , the detection rate of such sources are 
0 ^ 


expected to be roughly comparable to those end their jet 
activity at an early stage of their evolution (< 10 ® yr). 
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